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Abstract 



This paper concerns the reconstruction of an anisotropic conductivity tensor 7 from 
internal current densities of the form J = 7V1/, where u solves a second-order elliptic equation 
V • (7VU) = on a bounded domain X with prescribed boundary conditions. A minimum 
number of such functionals equal to n + 2, where n is the spatial dimension, is sufficient to 
guarantee a local reconstruction. We show that 7 can be uniquely reconstructed with a loss 
of one derivative compared to errors in the measurement of J. In the special case where 7 
is scalar, it can be reconstructed with no loss of derivatives. We provide a precise statement 
of what components may be reconstructed with a loss of zero or one derivatives. 

1 Introduction 

Hybrid medical imaging modalities are extensively studied in the bio-engineering community. 
Such methods aim to combine high-contrast, such as the one found in the modalities Electrical 
Impedance Tomography (EIT) or Optical Tomography (OT), with high-resolution, as is observed 
in the modalities Magnetic Resonance Imaging (MRI) or ultrasound. The high-contrast modality 
EIT aims to locate unhealthy tissues by reconstructing their electrical conductivity 7 from 
current boundary measurements. This leads to an inverse problem known as Calderon's problem. 
Extensive studies have been made on uniqueness properties and reconstruction methods for this 
inverse problem [33]. Unfortunately, the problem is severely ill-posed and yields images with 
poor resolution. 

It is sometimes possible to leverage a physical coupling between a high-contrast, low-resolution 
modality and a high-resolution, low-contrast modality. Such a coupling typically provides in- 
ternal functionals of the unknown coefficients of interest and greatly improve its resolution 
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[1, 2, 3, 6, 7, 22, 29, 32]. Different types of internal functionals, such as current densities 
and power densities, corresponding to different physical couplings have been analyzed to re- 
cover the unknown conductivity. In the case of power densities, we refer the reader to, e.g., 
[4, 5, 9, 16, 17, 22, 23, 24, 25]. 

In this paper, we consider the Current Density Impedance Imaging problem (CDII), also 
called Magnetic Resonance Electrical Impedance Tomography (MREIT) of reconstructing an 
anisotropic conductivity tensor in the second-order elliptic equation, 

n 

V • ( 7 V«) = £ ditf'dju) = (X), u\ dx =g, (1) 

from knowledge of internal current densities of the form H = "fVu, where u solves (1). To 
be consistent with earlier publications, where the notation H is used systematically to denote 
internal functionals, we use H to denote current densities rather than the more customary 
notation J. Here X is an open bounded domain with a C 2 ' a or smoother boundary dX. The 
above equation has real-valued coefficients and 7 is a symmetric tensor satisfying the uniform 
ellipticity condition 

K^Uf < £-t£ < k||£|| 2 , £eM n , for some k > 1, (2) 

so that (1) admits a unique solution in H l (X) for g G H^{dX). 

Internal current density functionals H can be obtained by the technique of current density 
imaging. The idea is to use Magnetic Resonance Imaging (MRI) to determine the magnetic 
field B induced by an input current /. The current density is then defined by H = V x B. We 
thus need to measure all components of B to calculate H , which may create some difficulties in 
practice, but this is the starting point of this paper. See [11, 30] for details. 

A perturbation method to reconstruct the unknown conductivity in the linearized case was 
presented in [12]. In dimension n = 2, a numerical reconstruction algorithm based on the 
construction of equipotential lines was given in [18]. Kwon et al [19] proposed a J-substitution 
algorithm, which is an iterative algorithm. Assuming knowledge of only the magnitude of 
only one current density \H\ = |7Vit|, the problem was studied in [26, 27, 28] (see the latter 
reference for a review) in the isotropic case and more recently in [10, 21] in the anisotropic case 
with anisotropy known. In [14, 20], Nachman et al. and Lee independently found a explicit 
reconstruction formula for visualizing log 7 at each point in a domain. The reconstruction with 
functionals of the form 7*Vit is shown in [15] in the isotropic case. For t = 0, the functionals 
are given by solutions of (1), then a more general complex- valued tensor in the anisotropic case 
was presented in [8]. In [31], assuming that the magnetic field B is measurable, Seo et al. gave 
a reconstruction for a complex-valued coefficient in the isotropic case. 

In the present work, we study the inverse problem in the anisotropic setting with a set of 
current densities Hj = jVuj for 1 < j < m, where Uj solves (1) with prescribed boundary 
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conditions gj. We propose sufficient conditions on tyi and the choice of {[Jj\ <.j<m 

such that the 

reconstruction of 7 is unique and satisfies elliptic stability estimates. 

2 Statement of the main results 

For X C M n , we denote by the set of conductivity tensors with bounded components 

satisfying the uniform ellipticity condition (2). Then for k > 1 an integer and < a < 1, we 
denote 

c£ a (X) := {7 G £(X)| lpq eC k > a (X), l<p<q<n}. 

In what follows, by "solution of (1)" we may refer to the solution itself or the boundary condition 
that generates it, i.e. g = u\qx £ H^(dX). We will consider collections of measurements of the 
form 

Ri : 7 Hi{i) = -fVui, 1 < i < m, (3) 

where i*j solves (1) with boundary condition g^. We decompose 7 into the product of a scalar 
factor /3 with an anisotropic structure 7 

7 := £7, /3 = (det7)«, det7 = 1. (4) 

Since 7 satisfies the uniform elliptic condition (2), (3 is bounded away from zero. 

From knowledge of a sufficiently large number of current densities, the reconstruction formu- 
las for (3 and 7 can be locally established in terms of the current densities and their derivatives 
up to first order. 

2.1 Main hypotheses 

We begin with the main hypotheses that allow us to setup a few reconstruction procedures. 

The first hypothesis aims at making the scalar factor /3 in (4) locally reconstructible via a 
gradient equation. 

Hypothesis 2.1. There exist two solutions (141,1*2) of (1) and Xq C X convex satisfying 

inf ^1(1*1,1*2) > Co > where T\ (1*1,1*2) := |Vi*i| 2 |Vi*2| 2 — (Vi*i • VU2) 2 . (5) 

xex 

On to the hypotheses for local reconstructibility of 7, we first need to have, locally, a basis 
of gradients of solutions of (1). 

Hypothesis 2.2. There exist n solutions (1*1, . . . ,u n ) of (1) and Xq C X satisfying 

inf JFiiui, . . . ,u n ) > Co > 0, where Jj^tti, • • • , u n ) := det(Vi*i, . . . , V« n ). (6) 

x£X 
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Let us now pick tti, ■ ■ ■ , u n satisfying Hyp. 2.2 and consider additional solutions {u n+ k} 
Each additional solution decomposes in the basis (Vui, . . . , Vu n ) as 



m 

k=V 



Vu n+fc = 22fi l k Vuj, 1 < k < m, (7) 

i=l 

where, as shown in [5] for instance, the coefficients [i\ take the expression 

i i 

det(Vui, . . . ,Vu n+b . . . ,Vu n ) det(i?i, . . .,H n+k , . . .,H n ) 



det(V«i, . . . , Vit n ) det(#i, . . . ,H n ) 



in particular, these coefficients are accessible from current densities. The subsequent algorithms 
will make extensive use of the matrix- valued quantities 

Zk = [Zk,i\ ■ ■ ■ \Z k , n ] , where Z k>i := Vfi l k , l<k<m (8) 

In particular, the next hypothesis, formulating a sufficient condition for local reconstructibility 
of the anisotropic part of 7 is that, locally, a certain number of matrices Z k (at least two) satisfies 
some rank maximality condition. 

Hypothesis 2.3. Assume that Hypothesis 2.2 holds for some (u±, . . . ,u n ) over Xq C X and 
denote by H the matrix with columns H\, . . . , H n . Then there exist u n +i, . . . , it n +m solutions of 
(1) and some X' C Xq such that the x -dependent space 

W := span{(Z k H T n) sym , Q G A n (R), 1 < k < m) C S n (R) (9) 

has codimension one in S n (M.) throughout X' . 

An alternate approach to reconstruct 7 is to set up a coupled system for u± , . . . , u n satisfying 
Hyp. 2.2 globally. This system of PDEs can be derived under the following hypothesis (part A). 
From this system and under an additional hypothesis (part B), we can derive an elliptic system 
from which to reconstruct u±, . . . ,u n . 

Hypothesis 2.4. A. Suppose that Hypothesis 2.2 is satisfied over Xq = X for some solutions 
(u±, . . . ,u n ). There exists an additional solution ii n +i of (1) whose matrix Z\ defined by 
(8) is uniformly invertible over X , i.e. 



inf det Z x > c > 0, (10) 



for some positive constant cq. 
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B. There exist n + 2 solutions ui, ... , it n +2 suc/i i/iai . . . ,u n ,u n+ 2) satisfy (A), and two 
A n (M) -valued functions £l\(x), ^(x) such that the matrix 

S = (Z^Z^U^x) + HZfn 2 {x)) sym {with Z% := Z 2 " T ) (11) 

satisfies the ellipticity condition (2). 

The first important result to note is that the hypotheses stated above remain satisfied under 
some perturbations of the boundary conditions or the conductivity tensor for smooth enough 
topologies. 

Proposition 2.5. Assume that Hypothesis 2.1, 2.2, 2.3 or 2. 4 holds over some Xq C X for a 
given number m of solutions of (1) with boundary conditions gi, . . . , g m . Then for any < a < 1, 
there exists a neighborhood of (gi, . . . ,g m -,l) open for the C 2,a (dX) m x C 1,a (X) topology where 
the same hypothesis holds over Xq. In the case of 2.4-B, it still holds with the same A n (JH) -valued 
functions Q\ and J^- 



2.2 Reconstruction algorithms and their properties 

Reconstruction of /3 knowing 7. Under knowledge of 7 and using two measurements H\,H2 
coming from two solutions satisfying Hyp. 2.1 over some Xq C X, we can derive the following 
gradient equation for log f3 

V1 °g/ ? = (l^il 2 d(7 _1 #i) - (Hi • H 2 ) d(T 1 H 2 )) (7^1, Wr'^i 

' 1 (12) 

-nri2 d(7 -1 fli)(7#i,0. x^ x o, 

I -nil 

where D := \Hi\ 2 \H2\' 2 — (Hi ■ H2) 2 is bounded away from zero over Xq thanks to Hyp. 2.1, and 
where the exterior calculus notations used here are recalled in Appendix A. 

Equation (12) allows us to reconstruct (3 under the knowledge of /3(xo) at one fixed point 
in Xq by integrating (12) over any curve starting from some xq S Xq. This leads to a unique 
and stable reconstruction with no loss of derivatives, as formulated in the following proposition. 
This generalizes the result in [14] to an anisotropic tensor. 

Proposition 2.6 (Local uniqueness and stability for j3). Consider two tensors 7 = (3^ and 
7' = (3'j', where 7,7' £ W 1,oc (X) are known. Suppose that Hypothesis 2.1 holds over the 
same Xq C X for two pairs (ui,U2) and (u^u^), solutions of (1) with conductivity 7 and 7', 
respectively. Then the following stability estimate holds for any p > 1 

|| log/3 - log/S'Hjyp.oop^) < e + C I ^2 \\ H i - #il|wi>.°°(X) + 117 - l'Wwp^(x) I (13) 

V=l,2 / 

Where eo = | log/?(xo) — log/3'(xo)| is the error committed at some fixed xq G Xq. 
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Algebraic, local reconstruction of 7: On to the local reconstruction of the anisotropic 
structure, we start from n -\- m solutions (iti, . . . , ii n +m) satisfying hypotheses 2.2 and 2.3 over 
some Xq C X. In particular, the linear space W C «S , n (M) defined in (9) is of codimension 
one in S* n (]R). We will see that the tensor 7 must be orthogonal to W for the inner product 
{A, B) := AijBij = tr (AB T ). Together with the conditions that det7 = 1 and 7 is positive, 
the space W, known from the measurements H\, . . . ,H n+m completely determines 7 over Xq. 
In light of these observations, a constructive reconstruction algorithm based on a generalization 
of the cross-product is proposed in section 4.2. This approach was recently used in [23] in the 
context of inverse conductivity from power densities. This algorithm leads to a unique and stable 
reconstruction in the sense of the following proposition. 

Proposition 2.7 (Local uniqueness and stability for 7). Consider two uniformly elliptic tensors 
7 and 7'. Suppose that Hypotheses 2.2 and 2.3 hold over the same Xq C X for two n + m-tuples 
{ui}"j~ 1 m and {^}" = t 1 m , solutions of (1) with conductivity 7 and 7', respectively. Then the 
following stability estimate holds for any integer p > 

||7 - l'\\wp<°°(x ) - ^ W Hi ~~ H 'i\\wp+ 1 -°°{x)- (14) 

i=l 

Joint reconstruction of (7,/?), stability improvement for V x 7 -1 . Judging by the 
stability estimates (14) and (13), reconstructing /3 after having reconstructed 7 is less stable 
(with respect to current densities) than when knowing 7. This is because in the former case, 
errors on W p,co -norm in 7 are controlled by errors in W p+1,co norm in current densities. In 
particular, on the W p,oc scale, stability on (3 is no better than that of 7, and joint reconstruction 
of (7, (5) using the preceding two algorithms displays the following stability, with 7 = ^7 

n+m 

||7 - l'\\wP.™{X ) < C ^ W H i ~ H' i \\wP+ 1 '°°(X)- ( 15 ) 
i=l 

However, once 7 is reconstructed, some linear combinations of first-order partials of 7" 1 can be 
reconstructed with better stability. These are the exterior derivatives of the columns of 7 _1 , a 
collection of n 2 (n — l)/2 scalar functions which we denote V x 7 -1 and is reconstructed via the 
formula 

8 ql pl -8 pl ql = H il ( 1 '"d p H ji -^dqHji), l<l<n, l<p<q<n, (16) 

derived in Sec. 4.3 and assuming that we are working with a basis of solutions satisfying 
Hypothesis 2.2. The stability statement (15) is thus somewhat improved into a statement of the 
form 

n+m 

||7 - l'\\wp>°°(Xo) + l|V x (7 _1 - 7 ~ 1 )|| W / P ,cx.( Xo ) < C ^2 W H i ~ H'iWwv+^ix), (I 7 ) 

i=l 
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where we have defined 



n 



V x (7 1 - 7' ^Hwrp.-^o) := J] S ~ ^7 j7 |kp.«(X )- 



J=l l<i<j<n 



Global reconstruction of 7 via a coupled elliptic system. While the preceding approach 
required a certain number of additional solutions, we now show how one can setup an alternate 
reconstruction procedure with only m = 2 additional solutions satisfying Hyp. 2.4. A microlocal 
study of linearized current densities functionals shows that this is the minimum number of 
functionals necessary to reconstruct all of 7. 

The present approach consists is eliminating 7 from the equations and writing an elliptic 
system of equations for the solutions u,-; see [5, 22, 23] for similar approaches in the setting of 
power density functionals. The method goes as follows. Assume that Hypothesis 2.2 holds for 
some (tii, • • • , u n ) over Xq = X and denote [VET] = [Vtii, • • • , Vu n ] as well as H = [Hi, • • • , H n ]. 
Since H = 7 [VET], we can thus reconstruct 7 by 7 = [VET] _1 .ff once [VET] is known. We now show 
that we may reconstruct [VET] by solving a second-order elliptic system of partial differential 
equations. 

When Hyp. 2. 4. A is satisfied for some ti n +i and considering an additional solution ii n +2 
and its corresponding current density, we first derive a system of coupled partial differential 
equations for (tii, . . . , u n ), whose coefficients only depend on measured quantities. 

Proposition 2.8. Suppose n + 2 solutions (u±, . . . , u n+ 2) satisfy Hypotheses 2.2 and 2.4-A and 
consider their corresponding measurements Hi = {flf}^ • Then the solutions (tii,-- - ,u n ) 
satisfy the coupled system of PDE 's 



for 1 < j < n and 1 < p < q < n, and where the vector fields vfj } only depend on the 
current densities Hi. 

If additionally, u n+ 2 is such that Hyp. 2.4.B is satisfied, we can deduce a strongly coupled 
elliptic system for (tii, ■ ■ ■ ■, u n ) from (18). 

Theorem 2.9. With the hypotheses of Proposition 2.8, assume further that Hypothesis 2.4-B 
holds for some A n (M) -valued functions 



Z%Z T (e p (g>e q -e q ® e p ) : V 2 uj + «?? • Vu; = 0, 
HZl{e p <g> e q - e q <g> e p ) : V 2 Uj + tf/ • Vu; = 0, 



Ujldx = 9j 



(18) 




l<p<<j<n 



Then (tii, ■ ■ ■ , 



u n ) can be reconstructed via the strongly coupled elliptic system 



V • (SVuj) + Wij -Vtij = 0, Uj\ax=9j, l<j<n, 



(19) 
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(20) 



l<p<5<n 



Moreover, if system (19) with trivial boundary conditions has only the trivial solution, u±, . . . ,u n 
are uniquely reconstructed. Subsequently, 7 reconstructed as 7 = HlVU]" 1 satisfies the stability 
estimate 



for data sets Hj,Hj close enough in H -norm. 
2.3 What tensors are reconstructible ? 

We now conclude with a discussion regarding what tensors are reconstructible from current 
densities, based on the extent to which Hypotheses 2.1-2.4 can be fulfilled, so that the above 
reconstruction algorithms can be implemented. 

Test cases. 

Proposition 2.10. For any smooth domain X C W 1 and considering a constant conductivity 
tensor 70, there exists a non-empty C 2 ' a -open subset of [H 2 (<9X)] n+2 of boundary conditions 
fulfilling Hypotheses 2.1-2.4 throughout X. 

The second test case regards isotropic smooth tensors of the form 7 = (3I n , where we show 
that the scalar coefficient f3 can be reconstructed globally by using the real and imaginary parts 
of the same complex geometrical optics (CGO) solution. The use of CGOs for fulfilling internal 
conditions was previously used in [4, 8, 25]. 

Proposition 2.11. For an isotropic tensor 7 = /3I n with f3 6 H 2 +3+e (X) for some e > 0, there 
exists a non-empty C 2 ' a -open subset of [Hz(dX)] 2 fulfilling Hypothesis 2.1 thoughout X. 

Thanks to Proposition 2.5, we can also formulate the following without proof. 

Corollary 2.12. Suppose 7 is a tensor as in either Proposition 2.10 or 2.11. Then, for any 
< a < 1, there exists a C ,a -neighborhood of '7 for which the conclusion of the same proposition 
remains valid. 



7 - j'\\l*(x) + ||V x (7- 1 - 7 -^l^x) < C\\Hi - H'j 



(21) 



8 



Push-forwards by diffeomorphisms Recall that for ^ : X — > VP(X) a V7 1,2 -diffeomorphism 
and 7 G S(X), we define \T/* 7 the conductivity tensor push-forwarded by ^ from 7 defined over 
by 

^ 7 := (|J 1I ,|- 1 D^. 7 .D^)o^- 1 , := detDI'. (22) 

We now show that, whenever a tensor is being push- forwarded from another by a diffeormor- 
phism, then the local or global reconstructibility of one is equivalent to that of the other, in the 
sense of the Proposition below. While the existence of in Y,(*$>(X)) merely requires that VP 
be a W^^-diffeomorphism, our results below will require that VP be smoother and that it satisfies 
the following uniform condition over X 

C^ 1 <\.h\< C* for some > 1. (23) 

Proposition 2.13. Assume that Hypothesis 2.1, 2.2, 2.3 or 2.4 holds over some Xq C X for 
a given number m of solutions of (1) with boundary conditions gi, ■ ■ ■ ,g m - For \t : X — > ^(X) 
a smooth diffeomorphism satisfying (23), the same hypothesis holds true over ^/(Xq) for the 
conductivity tensor ^7 with boundary conditions {g\ o vf -1 , . . . , g m o vt" 1 ). In the case of Hyp. 
2.4-B, it holds with the following A n (M.) -valued functions defined over *$>(X): 

*5>±fli := [D^ ■ Qi ■ D^] o \S>~ 1 and f*0 2 := [| J^\D^ ■ ft 2 • D¥} o tf -1 . (24) 

In contrast to inverse conductivity problems from boundary data, where the diffeomorphisms 
above are a well-known obstruction to injectivity, Proposition 2.13 precisely states the opposite: 
if a given tensor 7 is reconstructible in some sense, then so is \l/*7, and the boundary conditions 
making the inversion valid are explicitely given in terms of the ones that allow to reconstruct 7. 

Corollary 2.14. Suppose 7 is a tensor as in either Proposition 2.10 or 2.11 and ^> : X — > ^(X) 
is a diffeomorphism satisfying (23). Then the conclusion of the same proposition holds for the 
tensor \l/*7 over *&(X) and boundary conditions defined over d(^(X)). 

Generic reconstructibility. We finally state that any C 1,a smooth tensor is, in principle, 
reconstructible from current densities in the sense of the following proposition. This result uses 
the Runge approximation property, a property equivalent to the unique continuation principle, 
valid for Lipschitz-continuous tensors. 

Proposition 2.15. Let X C R n a C 2,a domain and 7 G C^ a (X). Then for any xq G X, there 
exists a neighborhood Xq C X of xq and n + 2 solutions of (1) fulfilling hypotheses 2.2 and 2.3 
over Xq . 
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Outline: The rest of the paper is structured as follows. Section 3 covers the preliminaries, 
including the proof of Proposition 2.5. Section 4 presents the derivations of the local recon- 
struction algorithms: Sec. 4.1 covers the local reconstruction of f3 and proves Proposition 2.6; 
Sec. 4.2 covers the local reconstruction of 7 and the proof of Proposition 2.7; Sec. 4.3 justifies 
equation (16); Sec. 4.4 discusses the global reconstruction of 7 via an elliptic system, with a 
proof of Propositions 2.8 and 2.9. Finally, Section 5 discusses the question of reconstructibility 
from current densities, with the proofs of Propositions 2.10, 2.11, 2.13 and 2.15. 

3 Preliminaries 

In this section, we briefly recall elliptic regularity results, the mapping properties of the current 
density operator and we conclude with the proof of Proposition 2.5. 

Properties of the forward mapping. In the following, we will make use of the following 
result, based on Schauder estimates for elliptic equations. It is for instance stated in [13]. 

Proposition 3.1. For k > 2 an integer and < a < 1, if X is a C k+1,a -smooth domain, then 
the mapping (5,7) >- > u, solution of (1) , is continuous in the functional setting 

C k ' a {dX) x Cj" 1,a (X) -> C k ' a (X). 

As a consequence, we can claim that, with the same k, a as above, the current density 
operator (</, 7) 1— > "fVu is continuous in the functional setting 

C k ' a (dX) x c£" 1,a P0 -> C fe ~ 1 ' Q (X). 
Moreover, this fact allows us to prove Proposition 2.5. 

Proof of Proposition 2.5. Fixing some domain Xq C X and using Proposition 3.1, it is clear 
that the mappings 

h : (C 2 ' a {dX)) 2 x C^ a (X) 3 (01,02,7) ^ mfJi(ui,«2), 

■AO 

h : (C 2 ' a (dX)) n x C^(X) 3 ( 9l , . . . , 5n , 7 ) M- MT 2 (UU ■ ■ ■ ,«n), 

with T\, JF2 defined in (5), (6), are continuous, so / 1 _1 ((0,oo)) and / 2 _1 ((0, 00)) are open, which 
takes care of Hypotheses 2.1 and 2.2. Further, Hypothesis 2.3 is fulfilled if and only if condition 31 
holds. Again, using Prop. 3.1, the mapping ^3 := infx B with B defined in (31) is a continuous 
function of (g u . . .,g n+m ,j) G (C 2 > a (dX)) n+m x C]f{X) so that /^((O.oo)) is open. 



10 



Along the same lines, Hypothesis 2. 4. A is stable under such perturbations because the map- 
ping 

(C 2 > a (dX)) n+1 x C^ a (X) 3 ( gi , . . . )5 n+i,7) ^ inf detZi, 

is continuous whenever u\,. . . ,u n satisfy (6) over X. Finally, fixing two A Tl (K)-valued functions 
Cli(x) and ^(x), Hypothesis 2.4.B is fulfilled whenever 

(91, ■ ■ ■ ,9n+2, 7) G p| ^((O.oo)), (25) 

i=l 

where we have defined the functionals, for 1 < i < n 

Sl : (C 2 > a (dX)) n+2 x B ( 91 , . . . , 5n+2 , 7 ) H- inf det{5 P3 }i< Pj3 <i, 

with S = {Sp t q}i< Pj q<.n defined as in (11). Such functionals are, again, continuous, in particular 
the set in the right-hand side of (25) is open. This concludes the proof. □ 

4 Reconstruction approaches 
4.1 Local reconstruction of (3 

In this section, we assume that 7 is known and with W l, °° components. Assuming Hypothesis 
2.1 is fulfilled for two solutions ui,u 2 over an open set Xq C X, we now prove equation (12). 

Proof of equation (12). Rewriting (3) as -j7 _1 -ffj = Vuj and applying the operator d(-). Using 
identities (49) and (50), we arrive at the following equation for log/3: 

Vlog^A(7- 1 ^i) = rf(7" 1 ^) 5 3 = 1,2. (26) 
Let us first notice the following equality of vector fields 

V log p A (T l )Hi (jH lr ) = (Vlog/3-7^i)(7 _1 ^i) - |#i| 2 Vlog/3, 

so that 



V log p = -i- (V log p ■ jHi ft" 1 ^ - — ^ V log p Air 1 (7#! , •) 



(Vlog/3 • 7-ffi)7 -1 iJi - T^dtf 1 ^)^, •)• 



|J^i| 2 V |#i| 2 
It remains thus to prove that 

(Vlog/3 -7#i) = 1 (|tfi| 2 d(7 _1 #i) - (H 1 -H 2 )d(T 1 H 2 )) (7^1,7^2) 



11 



which may be checked directly by computing, for j = 1,2 



d^H^Hu^Hz) = d log /3 A (r^'X^i, 7# 2 ) 

= (V log p ■ jH^Hj ■ H 2 - (V log /3 ■ jH 2 ) (Hj ■ H{). 

Taking the appropriate weighted sum of the above equations allows to extract (Vlog/3 • 7 -Hi), 
and hence (12). □ 



Reconstruction procedures for /3, uniqueness and stability. Suppose equation (12) 
holds over some convex set Xq C X and fix xo £ Xo. Equation (12) is a gradient equation 
Vlog/3 = F with known right-hand side F. For any x £ Xo, one may thus construct f3{x) by 
integrating (12) over the segment [xo,x], leading to one possible formula 

P{x) = P(x ) exp (x - x ) ■ F((l - t)x + tx) dtj , x G X . (27) 

Proof of Proposition 2.6. Since det7 = 1, the entries of 7 -1 are polynomials of the entries of 7, 
so that the entries of the right-hand side of (12) are polynomials of the entries of H±, i?2>7 an d 
their derivatives, with bounded coefficients. It is thus straightforward to establish that 

||Vlog/3 - VIog/3'|| £oo(Xo) < C(\\H - H'\\ w x^ [x) + ||7 - i'\\ w i^ {X )) (28) 

for some constant C. Estimate (13) then follows from the fact that 

|| log/3 - log/3'|| Loc(Xo) < I log/3(x ) - log/3'(x )| + A(X)|| V log /3 - V log P'\\ L °o {Xo) , 

where A.(X) denotes the diameter of X. □ 

One could use another integration curve than the segment [xo,x] to compute fi(x). In order 
for this integration to not depend on the choice of curve, the right-hand side F of (12) should 
satisfy the integrability condition dF = 0, a condition on the measurements which characterizes 
partially the range of the measurement operator. 

When measurements are noisy, said right-hand side may no longer satisfy this requirement, 
in which case the solution to (12) no longer exists. One way to remedy this issue is to solve 
the normal equation to (12) over Xo (whose boundary can be made smooth) with, for instance, 
Neuman boundary conditions: 

-Alog/3 = -V-F (X ), d v \ogp\ dXo =F-v, 

where v denotes the outward unit normal to Xo. This approach salvages existence while pro- 
jecting the data onto the range of the measurement operator, with a stability estimate similar 
to (13) on the H s Sobolev scale instead of the W s '°° one. 
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4.2 Local reconstruction of 7 

We now turn to the local reconstruction algorithm of 7. In this case, the reconstruction is 
algebraic, i.e. no longer involves integration of a gradient equation. In the sequel, we work with 
n + m solutions of (1) denoted {«j}™j 1 m , whose current densities {Hi = 7Vuj}™j 1 m are assumed 
to be measured. 

Derivation of the space of linear constraints (9). Apply the operator d(7 _1 ) to the 
relation of linear dependence 



Since the 2-form vanishes, by applying two vector fields 7e p , 7e p , 1 < p < q < n, we obtain, 



Notice that the above equation means ( ; yZ k )piH q i = ( : jZ k ) q iHpi, which amounts to the fact that 
jZ k H T is symmetric. This means in particular that jZ k H T is orthogonal to A n (R), and for 
any Q S A n (M), we can rewrite this orthogonality condition as 



where the last part comes from the fact that 7 is itself symmetric. Each matrix Z k thus generates 
a subspace of iS n (M) of linear contraints for 7. Considering m additional solutions, we arrive at 
the space of constraints defined in (9). 

Algebraic inversion of 7 via cross-product. We now show how to reconstruct 7 explicitely 
at any point where the space W defined in (9) has codimension one. We define the generalized 
cross product as follows. Over an iV-dimensional space V with a basis (ei, • • • , ejv), we define 
the alternating N — 1-linear mapping Af : V^ -1 — > V as the formal vector- valued determinant 
below, to be expanded along the last row 




1 < i < n. 



Z kji Aj- 1 H i = 0, where Z k)i :=Vfi%, A: = 1,2,... 



H q iZk,i ■ je p — HpiZk i ■ je g 



= tr (jZ k H T n) = tr (j T z k H T n) = 7 : z k H T n = 7 : (z k H T nyy m 



(29) 



(Vi,e 1 ) 



Af(V 1} --- 



V N -!) : 



1 



(30) 



det(ei, • • • 



e^) 



(V N -. 1 ,e 1 ) 



e N 
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Af(Vi,--- , Vjv-i) is orthogonal to V\,- ■■ , Vjv-i- Moreover, M{V\,--- , V/v-i) vanishes if and 
only if (Vi, ■ ■ ■ , Vn-i) are linearly dependent. 

With this notion of cross-product in the case V = S n (M), we derive the following reconstruc- 
tion algorithm for 7. Adding m additional solutions, we find that W can be spanned by f)W := 
n ^" 2 matrices whose expressions are given in (9), picking for instance {ei®ej— &j®ei}i<i<:j<n 
as a basis for A n (M). The condition that W is of codimension one over Xq can be formulated 
as: 

inf B(x) > ci > 0, B:= V |detJV(J)|», (31) 

where a{ng — l,j)W) denotes the sets of increasing injections from [l,ns — 1] to [l,{fW], and 
where we have defined jV(J) = N(M Il , • • • ,M/ ng _ 1 ), where N is defined by (30) with V = S n (R). 
Then under condition (31), W is of rank ns — 1 in 5 n (R). 

Whenever (Mi, . . . , M ns _i) are picked in W, their cross-product must be proportional to 7. 
The constant of proportionality can be deduced, up to sign, from the condition det 7 = 1 so we 
arrive at ±| det N{M\, ■ ■ ■ , M ng _i)| «7 = N(M\, • • • , M ng _i). The sign ambiguity is removed 
by ensuring that 7 must be symmetric definite positive, in particular its first coefficient on the 
diagonal should be positive. As a conclusion, we obtain the relation 

|detAA(/)|^7 = sign(AAn(/))AA(/), / e a(n s - 1, flW). (32) 

This relation is nontrivial (and allows to reconstruct 7) only if (Mi, . . . , M ng _i) are linearly 
independent. When codim W = 1 but (JW > ns — 1, we do not know a priori which ns — 1 
-subfamily of W has maximal rank, so we sum over all possibilities. Equation (32) then becomes 

sign(Mi(/M(/)=£ 7 , (33) 

Iev(n s -WW) 

with B defined in (31). Since B > c\ > over Xq, 7 can be algebraically reconstructed on Xq 
by formula (33), where M is defined by (30) with V = 5 n (R). 

Uniqueness and stability. Formula (33) has no ambiguity provided condition (31), hence 
the uniqueness. Regarding stability, we briefly justify Proposition 2.7. 

Proof of Proposition 2. 7. In formula (33), the components of the cross-products N{I) are smooth 
(polynomial) functions of the components of the matrices Z^H, which in turn are smooth func- 
tions of the components of {i?i}!=i OT and their first derivatives, and where the only term appear- 
ing as denominator is det(i^i, . . . , H n ), which is bounded away from zero by virtue of Hypothesis 
2.2. Thus (14) holds for p = 0. That it holds for any p > 1 is obtained by taking partial deriva- 
tives of the reconstruction formula of order p and bounding accordingly. □ 
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4.3 Joint reconstruction of (7, 0) and stability improvement 

In this section, we justify equation (16), which allows to justify the stability claim (17). Start- 
ing from n solutions satisfying Hypothesis 2.2 over Xq C X and denote H = {Hij}fj =1 = 
[Hi \ . . . \H n ] as well as H pq := (H~ 1 ) pg . Applying the operator d(j •) to both sides of (3) 
yields d( , y~ 1 Hj) = d(Vuj) = due to (49). Rewritten in scalar components for 1 < j < n and 
1 < p < q < n 

= d q (Y l H tj ) - d p {^ l H l0 ) = (d ql pl - d^ 1 )^ + j^d.Htj - ^ l d p H ly 

Thus (16) is obtained after multiplying the last right-hand side by H Jt , summing over j and 
using the property that Y^j=i HijH^ 1 = 5a. 

4.4 Reconstruction of 7 via an elliptic system 

In this section, we will construct a second order system for (ui, ■ ■ ■ , u n ) with n+2 measurements, 
assuming Hypotheses 2.2 and 2. 4. A hold with Xq = X. For the proof below, we shall recall the 
definition of the Lie Bracket of two vector fields in the euclidean setting: 

[X, Y) := (X • V)Y - (Y ■ V)X = (X'd^ej - (Y^X^-. 

Proof of Proposition 2.8. As is shown by (29), jZ^H 7 is symmetric. Multiplying both sides by 
7 _1 and using j~ l H = X7U, we see that Zk[X7U] T is symmetric. More explicitly, we have 

Zk,pidqUi = Zk^qidpUi, k = 1,2, (34) 

or simply Zk[VU] T = [VC/jZj. Assume Hypothesis 2. 4. A holds with Z2 invertible so that 
(^2,1) " • • j ^2,n) form a basis in M n . We define its dual frame such that Z$ j ■ Z2/1 = Denote 
Zjj = [Z21, • • • , Z\ J and Z\ = Z^ T ■ Then the symmetry of Z 2 [Vf/] T reads, 

Z$ d ■ Vm = Z% ti ■ Vuj, l<i<j<n. (35) 

Pick v a scalar function, we have the following commutation relation: 

{X • V)(Y • V)v = (Y ■ V)(X • V)v + [X, Y] ■ Vv 

Rewrite Zi tP id q = Zi )P ie q • V and apply Z\- • V to both sides of (34), we have the following 
equation by the above relations in Lie Bracket, 

[Z$j, Z hpi e q ] • Vui + (Z 1>pi e q • V)(Z$j ■ V) Ui = [Z\ y Z 1:qi e p ] • Vm + (Z 1:qi e p • V)(Z 2 %- • V)ttj 

(36) 
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where Zk,%j = Z^ : <g> e,-. Plugging (35) to the above equation gives, 

(Zi,pie q ■ V)(Z^i ■ V)uj + [Zjj, Z ltPi e q ) ■ Vuj = (Z^ qi e p ■ V)(Z^i ■ V)uj + [Z$j, Z^ qi e p ] ■ Vui 
Looking at the principal part, the first term of the LHS reads 

{Z 1)Pi e q • V)(Z^ • V) Uj = (Z^Zfe p ® e q ) : V\ + (Z 1>pi e q ■ V)^ • Vu^. 
Therefore, (36) amounts to the following coupled system, 

Z^Zj (e p ® e q - e q ® e p ) : V 2 Uj + ufj • Vu; = 0, Uj\ dx = gj, I < p < q < n (37) 

where 

■"f/ := [{Z\, p ie q - Z 1M e p ) ■ V] + [^j, Zi iPi e g - Z lsi e p ] . (38) 

Notice that H = 7[V£/] implies that -ff~ T [Vt/] T is symmetric. Compared with equation (34), 
we can see that the same proof holds if we replace Z2 by H~ T . In this case, the dual frame of 
H~ T is simply H. So (37) and (38) hold by replacing Z^ by H and defining u?? accordingly. □ 

We now suppose that Hypothesis 2.4.B is satisfied and proceed to the proof of Theorem 2.9. 

Proof. Starting from Hypothesis 2.4.B with A n (R)-valued functions of the form 

Sli(x) = ^2 uj l pq (x){e p ®e q -e q (Z)e q ), i = l,2, 

l<p<g<n 

we take the weighted sum of equations (18) with weights q , co p „. The principal part becomes 
S : V 2 Ui, which upon rewritting it as V • (SVuj) — (V • S) • Vui yields system (19). 

On to the proof of stability, pick another set of data H\ := {ff-}™^ 2 close enough to Hi in 
^1,00 norm) anc [ wr ite the corresponding system for u^, . . . , u' n 

-V • S'Vu'j + W[j ■ Vu'i = 0, 1 < j < n, (39) 

where <S" and W[a are defined by replacing Hi in (20) by Hj. Subtracting (39) from (19), we 
have the following coupled elliptic system for Vj = uj — u'j : 

-V • SVvj + Wij ■ Vvi = V ■ (S - S'^u'j + (Wlj - Wij) ■ VuJ, Vj \ dx = (40) 

The proof is now a consequence of the Fredholm alternative (as in [5, Theorem 2.9]). We recast 
(40) as an integral equation. Denote the operator Lq = — V • (5V) and define Lq 1 : H~ 1 (X) 3 
f 1— > v G Hq(X), where v is the unique solution to the equation 

-V-(SVv) = f (X), v\ dx =0. 
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By the Lax-Milgram theorem, we have IMI^ipq — CII/HiT-ipO? where C only depends on 
X and S. Thus Lq 1 : H^ 1 (X) — > H^(X) is continuous, and by Rellich imbedding, Lq 1 : 
L 2 {X) — > Hq(X) is compact. Define the vector space % = (i7g(JT)) n , v = (v\, . . . ,v n ), h = 
(Lq Vi, • • • , ioVn), where fj = V- (S-S^Vu'j + (W?- - Wij) • VttJ, and the operator V :U^U 
by 

P:H9v^Pv:= (Lq ^Wa • V^), • • • , i„>m ■ Vi>0) € ft. 

Since the are bounded, the differential operators Wjj-V : L^ — > L 2 are continuous. Together 
with the fact that Lq 1 : I? — > Hq is compact, we get that P : H. — > H. is compact. After applying 
the operator Lq 1 to (19), the elliptic system is reduced to the following Fredholm equation: 

(I + P)v = h. 

By the Fredholm alternative, if —1 is not an eigenvalue of P, then I+P is invertible and bounded 
||v|| w < IKl + Prl^Hhll^. Since L^ 1 : H-\X) H^{X) is continuous, h in (H%(X)) n is 
bounded by f =(/!,•••, / n ) in {H~ l (X)) n . 

\\M\U < W^WciH-^H^WH-^X)- 

Then we have the estimate, 

\\v\\u < ||(H-P) _1 ||£( W )||Lq ||£(i J --i ) H 1 )ll f lliJ-- 1 (X) 

Noting that Lq 1 is continuous and the RHS of (40) is expressed by Hj — Hj and their derivatives 
up to second order, we have the stability estimate 

ll u ~ u 'll//ipo ^ c \\ H i ~ H 'i\\m(x) 

where C depends on Hj but can be chosen uniform for Hj and Hj sufficiently close. Then 7 is 
reconstructed by 7 = LT[VC/] _1 and V x 7 _1 by (16), with a stability of the form 

117 - i'Lhx) + ||V x (7- 1 - y- 1 ))!!^*) < C\\Hi - H'tW^xy 

□ 

5 What tensors are reconstructible ? 

5.1 Test cases 

Constant tensors. We first prove that Hypotheses 2.1-2.4 can be fulfilled with explicit con- 
structions in the case of constant coefficients. 
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Proof of Proposition 2.10. Hypotheses 2.2 is trivially satisfied throughout X by choosing the 
collection of solutions U{(x) = X{ for 1 < i < n, then Hypothesis 2.1 is fulfilled by picking any 
two distinct solutions of the above family. 
Fulfilling Hypothesis 2.3. Let us pick 

Ui{x) : = Xi, 1 < i < n, 

1 _I " _i " 

u n+ i(x) : = -x T 7 2 >Jtj(e 3 - <8> ej)7o 2 ^ X,*?' = °j *p + k if ? / 9; 

2 j=i i=i (41) 

n-l 



1 -1 V— V — i 

n n+2 (3;) : = tt^To 2 z2( e j ® e i+i + e i+i ® e j)7 2x - 
In particular, H = and = V 2 « n +i for i = 1, 2, do not depend on a; and admit the expression 



_L _£ _L _i_ 

Z ^ = % 2 z2 t ^ e i ® e i) 7 o 2 and ^2 = 7 2 zJ( e J' ® e J+ 1 + e J+ 1 ® e i)^o 2 • 
i=i j=i 

We will show that the (x-independent) space 

W = span{(ZiF T 0)^ m , (2 2 ff T fi) sjm ,fie ^n(M)} 

has codimension one in S n (M.) by showing that W -1- C M70, the other inclusion D being evident. 

Let A E S n (R) and suppose that j4 _L W, we aim to show that A is proportional to 70. 

_ 1 _i 

The symmetry of AZ\H T implies that X^j=i^j e j ® e j7o 2 ^7o 2 ^ s symmetric. Denote -B = 
_ 1 _i 

7 2 A^y 2 £ S n (K) , we deduce that 

= tjBji, for 1 < i, j < n. 

Since B is symmetric and ti 7^ if i 7^ j, the above equation gives that Bij = for i 7^ j, 
thus -B is a diagonal matrix, i.e. B = Y17=l ^a e i ® e «- The symmetry of AZ2H T implies that 



X^jLi ( e j ® e j+i + e i+i ® e j)7o 2 ^7o 2 i s symmetric, which means that 

Bu(ej ® e j+ i + ej+i ® ej)(ei ® et) = ^ Bu (e* <8> e*) (e,- <g) e 3 - + i + e,+i <g> &j 

l<i<n l<i<n 
l<j<n-l l<j<n-l 

Write the above equation explicitly, we get 

n— 1 n—l 

X! B i+i«J+i e i ® e i+i + B .» e j • 1 <V X! % e J ® e J +1 + ''V I-/ - - 1 ® e i 

j=i i=i 
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Which amounts to 

n-l 

}^Bj +1> j +1 - Bjj)(e j+ i - e j+1 e,) = 
j=l 

Notice that {ej+i — ey+i0ej}i<j<„_i are linearly independent in A n (R), so Bj+ij+i = Bjj 
for 1 < j < ?7, — l, i.e. B is proportional to the identity matrix. This means that A must be 
proportional to 70 and thus W" 1- C K70. Hypothesis 2.3 is fulfilled throughout X. 



Fulfilling Hypothesis 2.4 with 7 = I n . We split the proof according to dimension. 

Even case n = 2m. Suppose that 77 = 2m, pick m = Xi for 1 < i < n, ii n +i = YlT=i x 2i-i%2i 
and 7j n +2 = Yl^Li — — ~- Then simple calculations show that 

m m 

Z\ = ^(e 2i _i 0e 2i + e 2 j 0e 2 j_i) and Z 2 = ^(e 2i _i e 2i _i - e 2i e 2i ). 



t=i 



t=l 



We have det Zi = (-l) m / so 2.4. A is fulfilled. Let us choose 

m 

fii := ^(e 2p e 2p _i - e 2p _i e 2p ) and f2 2 = 0, 
p=i 

then direct calculations show that 5 = (Z 2 Z^£l\ + HZfQ2) sym = Im which is clearly 
uniformly elliptic, hence 2.4.B is fulfilled. 

Odd case n = 3. Pick u[ = X{ for 1 < % < 3, 173+1 = a?ix 2 + £ 2 :e3 and u' s+2 = 2T[ x i + + 
2^x|, where ti,t 2 ,t 3 are to be chosen. In this case, H' = I3, Z[ = 2(ei e 2 + e 2 63) 
and (Z' 2 Y = Ya=i f i e i ® e * (note that Z£ fulfills 2.4.A). Pick n[(x) = e 2 e x - e x e 2 , 
f2 2 (a;) = e 2 e3 — e3 e 2 , simply calculations show that, 



s> = ((z' 2 r(z[) T n[(x) + H'(z[fn' 2 (x)) sy 



t 



Sa±l 
u 2 



-t 2 - 1 
fe+i 1 



(42) 



(ti,t 2 ,t^) must be such that S' is positive definite and tr {Z' 2 ) = (because u 2 solves (1)). 
This entails the conditions 



ti>0, ii(t 2 + l)<0, -(t 2 + l) ti 



> and t\ 



t 2 t 3 

t2 + t 3 ' 



These conditions can be jointly satisfied for instance by picking t\ = 6, t 2 = —2 and t 3 = 3, 
thus Hypothesis 2.4.B is fulfilled in the case n = 3. 
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Odd case n = 2m + 3. When n = 2m + 3 for m > 0, we build solutions based on the previous 
two cases. Let us pick 



1 < i < n, 

+1 = %2i-lX2i + X2m+lX2m+2 + ^2m+2^2m+3 



Ui — X{ , 
m 



i=l 



'"n+2 



i=i 



( X 2i-1 



Ad 



1 



i" x 2m+2 + g x 



"2m+3- 



Then one can simply check that Zj is of the form 



z 3 


02mx3 


03 x2m 


3 - 



J = 1,2, 



where ZjjZ'- are constructed as in the case n = 2m/n = 3, respectively. Accordingly, let 
us construct Ox,2 by block using the previous two cases, 





02m X 3 


03x2m 





and the S matrix so obtained becomes 



5 = (z\z1vl x + Hzfn 2 ^j 



sym 



l2m 


02mx3 


03 x2m 


S' 



where S' is the definite positive matrix constructed in the case n = 3. Again, Hypothesis 
2.4.B is fulfilled. 

Fulfilling Hypothesis 2.4 with 7 constant. Let {ui}^ denote the harmonic polynomials 
constructed in any case above (i.e. n even or odd) with 7 = I re , and denote Z®, Z® , H°, £1® 
and S° = (Z2*Zf T Q,® + H°Z® T Q2) sym the corresponding matrices. Define here, for 1 < % < n, 
Ui(x) := Vi(x) and for i = n+1, n+2, Ui(x) = Vi{^f~^x), all solutions of (1) with constant 7. Then 
we have that Z% = 7~2Z?7~2 for i = 1,2 and H = 7. Upon defining fij := 72Q972 g A n (R) 
for i = 1,2, direct calculations show that 

7*vTn 1 zjryTn \sym _ „,i c<0, 1 



5 = fix + n 2 )' 



72 5V- 



Whenever is non-singular, so is Z\ and whenever So is symmetric definite positive, so is S. 
The proof is complete. □ 
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Isotropic tensors. As a second test case, we show that, based on the construction of complex 
geometrical optics (CGO) solutions, Hypothesis 2.1 can be satisfied globally for an isotropic 
tensor 7 = f3I n when (3 is smooth enough. CGO solutions find many applications in inverse 
conductivity/diffusion problems, and more recently in problems with internal functionals [4, 25, 
8]. As established in [7], when j3 £ H% +3+£ (X), one is able to construct a complex-valued 
solution of (1) of the form 

u p = ±=e<>- x {l + i> p ), (43) 

where p £ C n is a complex frequency satisfying p ■ p = 0, which is equivalent to taking p = 
p(k + ik^) for some unit orthogonal vectors k, k -1- and p = \p\/V2 > 0. The remainder i/j p 
satisfies an estimate of the form pip p = 0(1) in C l {X). The real and imaginary parts of Vu p 
are almost orthogonal, modulo an error term that is small (uniformly over X) when p is large. 
We use this property here to fulfill Hypothesis 2.1. 

Proof of Proposition 2.11. Pick two unit orthogonal vectors k and k^, and consider the CGO 
solution u p as in (43) with p = p(k + ik.-^) for some p > which will be chosen large enough 
later. Computing the gradient of u p , we arrive at 

Vu p = e p - x (p + ip p ), with (f p : = V^ p - VpVlogy^, 

with sup^|<£>p| < C independent of p. Splitting into real and imaginary parts, each of which is 
a real- valued solution of (1), we obtain the expression 

Vu „ = pe^_ ^ + cog( ^ k± x) _ (k± + gin( ^ k± ^ 

V4 = ^ (V + p~ V?) cos(pk X • x) + (k + p~ VJ) sin(pk ± • *)) , 
from which we compute directly that 

2 2pk.-x 

|Vu«| 2 |V^| 2 - (VuJ • Vn^) 2 = (1 + o^ 1 )). 

Therefore, for p large enough, the quantity in the left-hand side above remains bounded away 
from zero throughout X, and the proof is complete. □ 



5.2 Push- forward by diffeomorphism 

Let : X — > ^f(X) be a VF^-diffeomorphism where X has smooth boundary. Then for 
7 £ the push-forwarded tensor ^^7 defined in (22) belongs to H(^f(X)) and ^ pushes 

forward a solution u of (1) to a function « = u o vf^ 1 satisfying the conductivity equation 

-V y ■ («VyV„v) = (#(X)), w| a( * ( x)) =<?°^ 1 , 
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moreover ^ and ^\ox induce respective isomorphisms of H l (X) and H^(dX) onto H 1 (^(X)) 
and Hv(d(V(X))). 

Proof of Proposition 2.13. The hypotheses of interest all fomulate the linear independence of 
some functionals in some sense. We must see first how these functionals are push-forwarded 
via the diffeomorphism ^ . For 1 < i < m, we denote V{ := ^+Ui = jjj o as well as 
fy±Hi := [*S>±j]V y Vi where y denotes the variable in ^f(X). Direct use of the chain rule allows 
to establish the following properties, true for any x € X: 

X7 Ul (x) = [D^] T (x)V y Vi(^(x)), 

Hi{x) = 7 Vni(x) = iJ^iD^-^H^ix)), (44) 

Zi(x) = [Dy] T (x)y*Zi{V(x)), 

where we have defined ^>-kZi the matrix with columns 

3 

[^ Zi ] 3 = _y / et ( V ^ ; ■ ■ ■ : ■ , 1<j<n . 

1 JJ v det(V^i,...,V^ n ) 

Hypotheses 2.1 and 2.2. Since [Dty] is never singular over X, relations (44) show that 
for any 1 < k < n, the vectors fields (Vui, . . . , Vitjt) are linearly dependent at x if and only if 
the vectors fields (V y vi, . . . ,V y Vk) are linearly dependent at ^(x). The case k = 2 takes care 
of Hyp. 2.1 while the case k = n takes care of Hyp. 2.2. 

Hypothesis 2.3. If we denote 

¥*W(*(a:)) = span{(^Z fc (^tf) T fl)^ m , n G A n (R),l< k < m} , 

direct computations show that 

W(x) = [D^{x)f • ¥*W(tf(s)) • [DV(x)\, 

thus since D*$>(x) is non-singular, we have that dimW(x) = dim *5>±W(*S>(x)), so the statement 
of Proposition holds for Hyp. 2.3. 

Hypothesis 2.4. The transformation rules (44) show that Z\ is nonsingular at x iff ty+Zi 
is nonsingular at ^(x), so the statement of the proposition holds for Hyp. 2. 4. A. 

Second, for two A n (M)-valued functions f2i(x) and ^(x), and upon defining ^Oi, ^>^2 as 
in (24), as well as 

tf*5 := ([**Z 2 ]- T [**Z 1 ] T **« 1 + [^ J ff][^Z 1 ] T ^fi 2 ) s?/m , 

direct use of relations (44) yield the relation 

S{x) = [D^ix)}- 1 ■ #*£(*(x)) • [D^(x)]~ T , xeX, 

and since is uniformly non-singular, S is uniformly elliptic if and only if ty+S is, so the 
statement of the proposition holds for Hyp. 2.4.B. □ 
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5.3 Generic reconstructibility 

We now show that, in principle, any C 1,Q -smooth conductivity tensor is locally reconstructible 
from current densities. The proof relies on the Runge approximation for elliptic equations, 
which is equivalent to the unique continuation principle, valid for conductivity tensors with 
Lipschitz-continuous components. 

This scheme of proof was recently used in the context of other inverse problems with internal 
functionals [8, 23], and the interested reader is invited to find more detailed proofs there. 

Proof of Proposition 2.15. Let xq £ X and denote 70 := 7(^0). We first construct solutions of 
the constant-coefficient problem by picking the functions defined in (41) (call them v±, . . . , v n +2) 
and by defining, for 1 < i < n + 2, u®{x) := v i(x) — fj(xo). These solutions satisfy V-(7oViij) = 
everywhere and fulfill Hypotheses 2.2 and 2.3 globally. 

Second, from solutions {u^l^t 2 , we construct a second family of solutions {t^}™!" 2 via the 
following equation 

V-( 7 V<) = (B 3r ), u r i \ aB3r =u i , l<i<n + 2, (45) 

where B 3r is the ball centered at x$ and of radius 3r, r being tuned at the end. The maximum 
principle as well as interior regularity results for elliptic equations allow to deduce the fact that 

lim max Wu^ — u®\\ r 2( Ro \ =0. (46) 

r-X) l<i<n+2 " 1 1 " C ^ 3r) K ' 

Third, assuming that r has been fixed at this stage, the Runge approximation property 
allows to claim that for every e > and 1 < i < n + 2, there exists gf E H?{dX) such that 

IK - <IU 2 (s 3t .) < e i where u i solves (!) with u t\dx = 9i, (47) 
which, combined with interior elliptic estimates, yields the estimate 

C C 

IK - u i \\c 2 (s;) ^ ^2 IK - <IU-(B 2r ) < 

Since r is fixed at this stage, we deduce that 

lim max ||< - <|| C 2 (B } = 0. (48) 

£^0 l<i<n+2 y 1 

Completing the argument, we recall that Hypotheses 2.2 and 2.3 are characterized by contin- 
uous functionals (say fi and ^3) in the topology of C 2,a boundary conditions. While the first step 
established that fi > and f% > for the constant-coefficient solutions, limits (46) and (48) tell 
us that there exists a small r > 0, then a small e > such that maxi<j< n+ 2 \\uf — v,® ||c 2 (s r (a:o)) 
is so small that, by the continuity of /2 and f%, these functionals remain positive. Hypotheses 
2.2 and 2.3 are thus satisfied over B r by the family {tif}™^ 2 which is controlled by boundary 
conditions. The proof is complete. □ 
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Remark 5.1 (On generic global reconstructibility) . Let us mention that from the local recon- 
structibility statement above, one can establish a global reconstructibility one. Heuristically, by 
compactness of X, one can cover the domain with a finite number of either neighborhoods as 
above or subdomains diffeomorphic to a half-ball if the point xq is close to dX , over each of which 
7 is reconstructible. One can then patch together the local reconstructions using for instance a 
partition of unity, and obtain a globally reconstructed 7. The additional technicalities that this 
proof incurs may be found in [8j. 

As a conclusion, for any C 1,a -smooth tensor 7, there exists a finite N and non-empty open 
set O C (C 2 ' a (dX)) N such that any {gi}f = i £ O generates current densities that reconstruct 7 
uniquely and stably (in the sense of estimate (17)) throughout X. 

A Exterior calculus and notations 

Throughout this paper, we use the following convention regarding exterior calculus. Because we 
are in the Euclidean setting, we will avoid the flat operator notation by identifying vector fields 
with one- forms via the identification e« = e* where {ei}f =1 and {e*}™ =1 denote bases of M n and 
its dual, respectively. In this setting, if V = V l ei is a vector field, dV denotes the two-vector 
field 

dV = {diV j -djV^Aej. 

l<i<j<n 

A two-vector field can be paired with two other vector fields via the formula 

A A B(C, D) = {A- C){B ■ D) — (A ■ D)(B ■ C), 
which allows to make sense of expressions of the form 

dV(A, ■)= ( diVJ ~ d j yl )(( A ■ e *) e i " ( A ■ e J') e *)- 

l<i<j<ri 

Note also the following well-known identities for / a smooth function and V a smooth vector 
field, rewritten with the notation above: 

d(V/) = 0, feC 2 (X), (49) 
d(fV) = V/ A V + fdV. (50) 
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